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ABSTRACT 

We present sub-arcsecond resolution imaging of the A = 2.7 mm continuum emission 
from the young, embedded system L1551 IRS5 using the nine-element, high-resolution 
configuration of the BIMA array. The observed emission arises from two compact 
sources separated by 0Y35, coinciding with the two sources seen at A = 2 cm and 
A = 1.3 cm. When the high resolution data is combined with data from two compact 
configurations, L1551 IRS5 is argued to consist of a protobinary system separated by 
~50 AU with individual circumstellar disks, a circumbinary structure, and a large-scale 
envelope. The characteristic masses of the components are: 0.024 Mq for the northern 
circumstellar disk, 0.009 Mq for the southern circumstellar disk, 0.04 Mq for the 
circumbinary material, and 0.28 Mq for the envelope. 



Subject headings: starsxircumstellar — stars:formation — binaries: close — stars: 
individual (L1551 IRS5) — Infrared: stars — Radio Continuum: stars 



1. Introduction 



First detected in an infrared survey of the L1551 cloud (Strom, Strom, & Vrba 1976), 
L1551 IRS5 is a prototypical young stellar system, with a strong bipolar molecular outflow (Snell, 
Loren, & Plambeck 1980), an optical jet (Mundt & Fried 1983), HH objects (Herbig 1974), and 
an envelope-disk structure in the surrounding material (Keene & Masson 1990). Located at a 
distance of 140 pc (Elias 1978) and exhibiting a luminosity of ~ 28Lq (Butner et al. 1991), 
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L1551 IRS5 was one of the denning examples for Class I sources in the classification scheme of 
Adams, Lada, and Shu (1987) and has been used as an archetype in the current paradigm for 
single-star formation (Shu et al. 1993). But is it really a single-star system? 

High resolution A = 2 cm continuum observations of L1551 IRS5 show two compact sources 
with a separation of ~ 0'.'28 (Bieging & Cohen 1985; Rodriguez et al. 1986) which have been 
interpreted as either a protobinary system (Bieging & Cohen 1985), or the inner ionized edges 
of a gas and dust toroid surrounding a single star (Rodriguez et al. 1986). The latter is the 
most widely accepted interpretation, but comparisons with A = 2 cm emission from other young 
binary systems such as T Tau and Z CMa (Bieging, Cohen, & Schwartz 1984; Schwartz, Simon, & 
Zuckerman 1983), suggest that the binary interpretation is also viable. 

Under the assumption that L1551 IRS5 is a single star system, Keene and Masson (1990) 
modeled A = 2.7 mm interferometric observations to deduce the presence of a 45 AU radius 
circumstellar disk within an envelope. This envelope, which extends out ~1000 AU from IRS5, 
contains 0.1 to 1 Mq of material (Ladd et al. 1995; Fuller et al. 1995). High resolution 
JCMT-CSO interferometric observations at A = 870 fim resolved the compact central emission 
(Lay et al. 1994), and the emission was modeled as arising from an 80 AU radius Gaussian source, 
inferred to be an accretion disk around the young star. 

In this letter, we present sub-arcsecond imaging of the A = 2.7 mm continuum emission from 
the L1551 IRS5 system. These observations re-open questions about the binarity of the system 
and the distribution of the surrounding material. 

2. Observations and Data Reduction 

L1551 IRS5 was observed in three array configurations of the 9-element BIMA ArrayQ (Welch 
et al. 1996). The longest baselines were 1 km N-S and 900 m E-W, yielding a maximum projected 
baseline of 480 kA (1.4km); the shortest baselines were limited by the antenna size of 6.1 m, 
yielding a minimum projected baseline of 2.2 kA. This range in projected baselines provides images 
with a minimum resolution of 0'.'3, fully sampled to sizes as large as 60". 

For the high resolution configuration (March 1, 1996), atmospheric phase fluctuations were 
tracked by switching the antennas between source, phase calibrator, and a nearby weak quasar 
on a two minute cycle. The usefulness of this quick switching technique has been demonstrated 
at the VLA (Holdaway & Owen 1995). The main phase calibrator (0530+135) was used to track 
rapid atmospheric phase fluctuations. The secondary quasar (0449+113) was used to track slow 
phase drifts due to the difference in airmass between the primary calibrator and source and, more 
importantly for this array, phase drifts due to uncertainties in baseline length. 



4 The BIMA Array is operated by the Berkeley Illinois Maryland Association under funding from the National 
Science Foundation. 
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The digital correlator was configured with two 700 MHz bands centered at 107 GHz and 
109 GHz. The flux amplitude calibration assumed a flux of 6.8 Jy for 0530+135, as observed 
in the following month's compact array. The coherence of the atmosphere was checked on the 
quasars; the uncertainty in the amplitude calibration is 20%. Absolute positions in our map have 
uncertainty due to the uncertainty in the antenna baselines and the statistical uncertainty from 
the signal-to-noise of the observation. These two factors add in quadrature to give an absolute 
positional uncertainty of / /14. The lower resolution data (acquired on October 3, 1996, February 
2, 1997, and March 8, 1997) used 0530+135 to track phase variations and Mars for amplitude 
calibration. 

The L1551 IRS5 data were imaged in four ways which stress structures present on different 
spatial scales. Figure 1 shows four maps: two with robust weightings of the visibilities (robust = 0.5 
yielding a 3'.'25 x 3'.'04 beam and robust = -0.25 yielding a l'.'ll x / /84 beam), one with natural 
weighting of only the high resolution A array data restored with the fitted "clean" beam 
(0'.'73 x 0'.'31 beam), and one with the A array data restored with a circular tf!31 "clean" beam. 
The latter technique strongly emphasizes the high resolution information present in the A array 
u,v data. With maximum projected baselines ranging from 320 kA to 480 kA, the smallest fringe 
spacings in our dataset ranges from / . / 64 to 0'/44; hence information down to size scales of / /2 to 
Cf.'3 is present in the u,v data. High resolution maps of the secondary quasar 0449+113, using the 
standard technique and using the / /31 "clean" beam, were consistent with a point source. 

3. Results 

Figure la (3" resolution) has a peak flux of 122+3 mJy beam -1 , and the integrated flux in 
a 8" box centered on the source is 162+6 mJy. A Gaussian fit to the image gives a deconvolved 

o 

Gaussian source size of l'/78 x l'/75 and PA=68 . Figure lb (1" resolution) has a peak flux of 
78+3 mJy beam -1 , and the integrated flux in a 3" box centered on the source is 143+10 mJy. A 

o 

Gaussian fit to the image gives a deconvolved Gaussian source size of 0'.'92 x (y/61 and PA=157 . 
Figure lc shows the map of the A array data alone restored with the Gaussian fitted clean beam. 
The peak flux in the map is 45+5 mJy beam -1 , and the integrated flux in a l'/3 box centered 
on the source is 75+11 mJy. Although it is not obvious in Figure lc, over \ of the flux present 
in the lowest resolution map is now gone and the peak flux is roughly | of that in Figure la. 
Despite the elongated "clean" beam, the remaining emission is clearly extended north-south in 
the CLEANed image; a Gaussian fit to the image gives a deconvolved Gaussian source size of 

o 

C/53 x / .'32 and PA=7.2 . Figure Id shows the A array data restored with the circular / /31 beam. 
The north-south extension is obvious in this map and there is no hint of east-west extension. The 
peak flux is 38 mJy beam -1 corresponding to a brightness temperature of 41 K. 

The images in Figure la and lb emphasize the overall emission from the L1551 system. 
The reconstructions in Figure lc and Figure Id highlight the small scale emission which is more 
compact than expected for the disk size estimates of Keene and Masson (1990) and Lay et al. 
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(1994). The compact emission is consistent with arising from two point sources, as seen at 
A = 2 cm and A = 1.3 cm (Rodriguez et al. 1986; Koerner & Sargent 1997). A two Gaussian fit to 
the A = 2.7 mm emission in Figure Id yields the following positions hereafter labeled IRS5 A and 
IRS5 B: IRS5 A: a(J2000) = 04 fe 31 m 34M43, c5(J2000) = 18° 08 05'.'09 and IRS5 B: a(J2000) = 
04 A 31 m 34M41, <5(J2000) = 18° 08' 04'/ 74. These positions agree to within 0'/05 with the A = 1.3 cm 
source positions of Koerner and Sargent (1997). The separation of the two sources is (X'35, 
corresponding to 49 AU. Both sources have deconvolved sizes of < 0?3. A two point source fit 
yields flux densities of 45±6 mJy for IRS5 A and 23±6 mJy for IRS5 B. The total flux density in 
the compact sources is then 68±9 mJy. 

4. Comparisons with Centimeter High Resolution Data 

High resolution centimeter wavelength images of L1551 IRS5 show two point-like sources and 
an extended jet (Bieging & Cohen 1985). The jet is detected only at long centimeter wavelengths; 
the two point sources dominate the flux at shorter wavelengths. The A = 2 cm flux densities are 
1.2 mJy for IRS5 A and 0.93 mJy for IRS5 B (Rodriguez et al. 1986). Recent VLA observations 
also resolved the two sources at A = 1.3 cm (Koerner & Sargent 1997) and yielded flux densities 
of 2.0±0.2 mJy and 1.5±0.2 mJy, respectively. The spectral indices between A = 2.0 and 1.3 cm 
are then cxa ~1.25 and ctB ~1.04, consistent with a ~1 estimated by Bieging and Cohen (1985). 
Extrapolating to 109 GHz, this emission could contribute as much as ~14.4 mJy and ^7.8 mJy, 
respectively, to the observed fluxes. Hence, the A = 2.7 mm flux is dominated by dust emission. 

The proposal of Rodriguez et al. (1986) that the A = 2 cm emission traces the ionized inner 
edge of a larger dusty torus is not consistent with the observed compact A = 2.7 mm emission. 
Since the millimeter emission directly probes the dust, we should easily see the torus in our 
high resolution maps. If there were a torus, the A = 2.7 mm emission would extend beyond the 
A = 2 cm sources and, in fact, peak outside of them. The binary interpretation of Bieging & 
Cohen is consistent with our image if the A = 2.7 mm emission arises from circumstellar disks 
within the binary system, while the A = 2 cm emission traces ionized gas associated with stellar 
winds or jets. 

5. The Structure of the L1551 IRS5 System 

Combining all observations to date, the L1551 IRS5 system consists of three main circumstellar 
components: a large-scale envelope (Keene and Masson 1990; Ladd et al. 1995), a disk or extended 
structure with a size scale of ~1" (Lay et al. 1994; Keene and Masson 1990), and an inner binary 
system as argued in section 4. How do these components fit together? To answer this question, 
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we compare our u,v data binned in annuli with simulated observations of models for the system, 
binned similarly. In the following subsections we discuss each component and derive characteristic 
masses. 



5.1. Binary Circumstellar Disks 

Figure 2a compares our u,v data with the Gaussian model from Lay et al. (1994) scaled to 
match the A = 2.7 mm flux at 50 kA. Above 100 kA, the Gaussian model is resolved out and does 
not fit the data; below 20 kA the data diverge from the model due to flux from the envelope. 
Figure 2b shows a two point source model with the separation and amplitudes given in section 3. 
The two point sources beat together to cause the variations in flux seen past lOOkA. The proposed 
binary disk system is evident only in our data; its separation is too small to be resolved in the data 
of Lay et al. (1994) or Keene and Masson (1990). In fact, due to the small angular size and the 
embedded nature of the binary system, the properties of the proposed disks are poorly constrained 
by observations to date. The projected separation and extent of the A = 2.7 mm emission suggests 
a maximum outer radius of 25 AU for the disks. To estimate the masses of the disks, we assume a 
standard power-law disk with parameters characteristic of the HL Tau disk, T^isk = 330(^^-) ' 5 
and Srfjsfe oc r -1 (Mundy et al. 1996; Beckwith & Sargent 1991). For dust properties, we adopt 
k=0.1( 120^ GHz ) CU]2 S -1 ' wmcn is consistent with other recent works (e.g. Osterloh and Beckwith 
1995; Ohashi et al. 1991; Beckwith & Sargent 1991). With these assumptions, the disk masses are 
M A ~0.024 M and M B ~0.009 M . 



5.2. The Envelope 

The excess emission in our robust weighted maps (Figures la and lb) compared to our highest 
resolution map (Figure Id) and the rise in flux on baselines shorter than 15 kA (Figure 2), are 
due primarily to the extended envelope. Our flux densities in the larger beams are consistent with 
previous measurements at similar resolutions: Keene and Masson (1990) find a peak flux of 130 
mJy beam -1 at A = 2.73 mm in a 2'/6 beam and a total flux of 150 mJy; Ohashi et al. (1996) 
measure a total flux of 160 mJy at A = 2.73 mm using a 4'.'5 beam. Our u,v data in Figure 2 and 
Figure 2 of Keene and Masson (1990), show similar fluxes around 10 kA, but our data has 20% to 
30% less flux from 40 kA to 70 kA. These differences are within the calibration uncertainties. 

The differences in flux densities at different resolutions, or equivalently the drop in flux density 
with u,v distance, can be used to estimate the properties of the envelope. Our data are broadly 
consistent with the envelope parameters determined by Ladd et aZ.(1995) and Fuller et aZ.(1995). 
Fitting the drop in flux between 2.6 kA and 15 kA with a power-law envelope model (p(r) oc r~ L5 
and T{r) oc r -0,5 ) combined with the two point source model from section 5.1, reasonable results 
are obtained for an envelope mass of ~0.44 Mq, an outer radius of ~1300 AU, and an inner 
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envelope radius of 30 AU (Figure 2c). Steeper envelope density laws (p(r) oc r 2 ) also fit the data 
with a characteristic mass and outer radius of 0.43 Mq and 1800 AU, respectively. 

5.3. The Circumbinary Structure 

Finally, an intermediate-sized structure, perhaps a circumbinary disk such as seen around 
GG Tau (Dutrey, Guilloteau, Simon 1994) or a "pseudo-disk" (Galli and Shu 1993), is needed 
to account for the structure resolved by Lay et al. (1994) and the compact structure deduced by 
Keene and Masson (1990). In our data, this structure is evident as the excess emission between 
30 kA and 90 kA in Figure 2c. As shown in Figure 2d, this excess can be fitted with a Gaussian 

o 

model consistent with that of Lay et al. (172 x / /7 PA = 160 ) with a flux of 30 mJy plus an 
envelope model with a mass of 0.28 Mq (p(r) oc r -15 ) and a radius of 1100 AU. The parameters 
of the envelope and the circumbinary structure are interdependent and hence only crudely 
determined. If the circumbinary structure has dust properties similar to the envelope parameters 
in section 5.2, the circumbinary structure has a rough mass of 0.04 Mq. 

To test the consistency of the above model with the Lay et al. (1994) data, we fit two different 
source structures to their A = 870 pm data: a single elliptical Gaussian (a single circumstellar 
disk) and a single elliptical Gaussian with two central point sources (a circumbinary disk with two 
small circumstellar disks), following the fitting procedure of Lay et al. (1994; also see Lay 1994). 
The model did not include envelope emission since the JCMT-CSO baselines ranged from 50 kA 
to 200 kA, where the envelope emission is completely resolved out. The single elliptical Gaussian 
model fits the A = 870 fim data very well, with parameters comparable to those found by Lay et 
al. (1994). The addition of two point sources to the single Gaussian model produces as good a fit 
as the single Gaussian model, but the FWHM of the Gaussian increases slightly. Hence, the data 
cannot distinguish between the single Gaussian and single Gaussian with point source models. If 
the circumbinary material is optically thick at A = 870 ^m, the Lay et al. data would not even 
see the embedded circumstellar disks. If the circumbinary material is not optically thick, the Lay 
et al. data place a limit on the flux from the circumstellar disks: at a 95% confidence level the 
circumstellar disks emit < 1.3 Jy at A = 870 pm. 

6. Young Binary systems 

Our data present the first direct detection of a close, embedded binary system. Proposed 
wider binary systems have been identified among embedded sources, e.g. IRAS 16293-2422 
(Wootten 1989), NGC 1333 IRAS4 (Sandell et al. 1991; Lay et al. 1995), and L1527 (Fuller, Ladd, 
& Hodapp 1996), but the number of such systems is actually quite small compared to the number 
of known embedded sources. Surveys of pre-main sequence (PMS) stars find that binary systems 
are at least as common among young visible stars as among main-sequence stars (Simon et al. 
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1992; Ghez, Neugebauer & Matthews 1993; Leinert et al. 1993; Reipurth & Zinnecker 1993); so 
binaries should be common among young, deeply embedded systems. That they have not often 
been seen is probably due to the lack of sub-arcsecond resolution observations which are necessary 
to resolve close binaries. The separation of the L1551 IRS5 system is near the median separation 
for main sequence binaries (~30 AU, Duquennoy & Mayor 1991). The low detection rate of 
wide embedded binaries is in rough agreement with the fraction of main sequence binaries with 
separations between 300 and several 1000 AU. 

L1551 IRS5 also ranks as one of the few close binary systems with significant dust emission 
associated with both components. Submillimeter wavelength surveys have generally found lesser 
amounts of dust emission associated with PMS binary systems than with young single stars 
systems (Simon et al. 1992, 1995). In a statistical comparison of binaries and single stars, Jensen, 
Mathieu & Fuller (1996) found that binaries with separations < 50-100 AU statistically have lower 
submillimeter fluxes than wider binaries, but wide binaries are indistinguishable from single stars; 
hence, the L1551 IRS5 system may be unusual. However, these studies concentrate on T Tauri 
stars and exclude the youngest sources, Class I or younger. It is possible that embedded close 
binaries, which are still accreting mass, have substantial circumstellar or circumbinary disks which 
disappear later when the envelope is no longer feeding-in material. 

7. Conclusions 

Sub-arcsecond A = 2.7 mm observations of L1551 IRS5 have resolved a compact central 
structure, which is most plausibly interpreted as a young binary system. The A = 2.7 mm 
continuum emission shows two peaks which are similar, in absolute position and separation, to the 
free-free emission observed at centimeter wavelengths. Our interpretation is that we are detecting 
thermal dust emission from small disks around the individual stars in a binary system and that 
the centimeter emission arises in the associated stellar winds. We propose that the L1551 IRS5 
system is composed of two circumstellar disks, located inside a circumbinary structure, embedded 
in a large-scale envelope. Simple modeling yields masses for these components: circumstellar 
disk masses of 0.024 Mq and 0.009 Mq for the northern and southern sources respectively, a 
circumbinary structure mass of 0.04 Mq, and an envelope mass of 0.28 Mq. The binary separation 
for L1551 IRS5 is about 50 AU, close to the median separation for main sequence binaries. The 
small number of young embedded binaries detected to date, probably reflects the inadequate 
angular resolution available in the earlier studies, rather than an intrinsic sparsity of binaries. 

We thank the Hat Creek staff for their efforts in the construction and operation of the long 
baselines array. We also thank the referee, Jocelyn Keene, and Oliver Lay for their comments. 
This work was supported by NSF Grants NSF-FD93-20238 and AST-9314847. LGM acknowledges 
support from NASA grant NAGW-3066. 
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Fig. 1. — A = 2.7 mm maps of the continuum emission from L1551 IRS5. a) Robust weighting of 

o 

0.5 map made with data from three arrays. The beam is 3. 25 x 3. 04 PA = 29 , and the RMS noise 
is 2.5 mJy beam" 1 . The contours are -3,-2,2,3,4,5,6,7,8,9,10,15,20,25,30,35 times 3.3 mJy beam" 1 
(the RMS from Panel b). b) Robust weighting of -0.25 map made with data from three arrays. 

o 

The beam is l'.'ll x 0'.'84 PA = 60 , and the RMS noise is 3.3 mJy beam . The contours are the 
same as in Panel a. c) Naturally weighted map made from only the A array data. The beam is 
C/73 x / /31 PA=47 and the RMS is 4.5 mJy beam -1 . The contours are in steps of 1 a starting at 
±2 a. d) A array naturally weighted data, restored with a circular C//31 beam. The contours and 
RMS are the same as in Panel c. The two crosses in Panels c and d mark the A = 1.3 cm source 
positions from Koerner & Sargent 1997. The restoring beam in each panel is shown in the lower 
left-hand corner. 
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Fig. 2. — The measured A = 2.7 mm visibilities binned in annuli (open squares) compared with 
different model visibilities (gray, closed squares), a) The Lay et al. (1994) model Gaussian scaled 
to match the A = 2.7 mm fluxes around 50 kA. b) Two point source model constrained by fitting 
Figure Id. c) Characteristic model fit with an envelope (0.44 Mq, p(r) oc r _L5 , T(r) oc r~ ' 5 , and 
1300 AU radius) and the two point sources from Panel b. d) Characteristic fit for a model with an 
envelope (0.28 M , p(r) oc r" 1 - 5 , T(r) oc r" - 5 , and 1100 AU radius), a Gaussian (30 mJy, l'/2 x 0"7 

o 

PA = 160 ), and the two point sources from Panel b. 
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